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Nanoprecipitates form during nucleation of multiphase equilibria in phase segregating multicom- 
ponent systems. In spite of their ubiquity, their size-dependent physical chemistry, in particular 
at the boundary between phases with incompatible topologies, is still rather arcane. Here we use 
extensive atomistic simulations to map out the size-temperature phase diagram of Cu nanoprecipi- 
tates in a-Fe. The growing precipitates undergo martensitic transformations from the body-centered 
cubic (BCC) phase to multiply-twinned 9R structures. At high temperatures, the transitions ex- 
hibit strong first-order character and prominent hysteresis. Upon cooling the discontinuities become 
less pronounced and the transitions occur at ever smaller cluster sizes. Below 300 K the hysteresis 
vanishes while the transition remains discontinuous with a finite but diminishing latent heat. This 
unusual size-temperature phase diagram results from the entropy generated by the soft modes of 
the BCC-Cu phase, which are stabilized through confinement by the a-Fe lattice. 

PACS numbers: 64.70.kd, 64.70.Nd, 64.75.Jk, 81.30.Hd 



The Fe-Cu system is a prototypical structurally mis- 
matched system with very small solubility of Cu in Fe. 
At low temperatures Fe adopts the body-centered cubic 
(BCC) phase, while Cu prefers the face-centered cubic 
(FCC) phase with BCC-Cu being mechanically unstable. 
Phase segregation in Fe-Cu alloys also has been and con- 
tinues to be the subject of considerable technological in- 
terest since Cu impurities contribute to high-temperature 
cmbrittlement of ferritic steels [3]. Based on a series of 
careful experimental studies three stages in the transfor- 
mation of Cu precipitates have been identified (TJ [51 0HH] : 
(i) Cu precipitates nucleate in the BCC structure of the 
host lattice, (ii) as they grow they undergo a marten- 
sitic phase transition to a multiply-twinned 9R structure, 
which (Hi) eventually transforms into FCC. In this Letter 
using a novel combination of molecular dynamics (MD) 
and Monte Carlo (MC) simulations as well as empirical 
interatomic potentials, we develop a comprehensive pic- 
ture of the correlation between precipitate structure and 
size as a function of temperature for the first two stages 
of the transformation process. 

The Fe-Cu system was modeled using semi-empirical 
interatomic potentials [9TI11] that provide an accurate de- 
scription of the alloy phase diagram and yield very good 
agreement with density-functional theory calculations for 
the energies of small Cu clusters in BCC-Fe. Simulations 
were carried out using the hybrid MC/MD algorithm in- 
troduced in [121 113) implemented in the massively par- 
allel MD code lammps [14] . Orthorhombic simulation 
cells were employed with cell vectors oriented along [111], 
[110], and [112]. Most simulations used simulation cells 
containing 42 x 52 x 45 unit cells (786,240 atoms) with 
additional simulations of smaller precipitates at low tem- 
peratures using 31 x 38 x 33 unit cells (310,992 atoms). 
Periodic boundary conditions were applied in all direc- 
tions, and pressure and temperature were controlled us- 



ing Nose-Hoover thermostat and barostat [15] . Simula- 
tions were run for at least 1.6 x 10 6 and up to 4 x 10 6 
MD steps using a timestep of 2.5 fs. Every 20 steps the 
MD simulation was interrupted to carry out a number 
of VC-SGC MC trial moves corresponding to 10% of a 
full MC sweep. Each simulation run thus comprises be- 
tween 8,000 and 20,000 attempts per particle to swap the 
atom type. The atomic structures were analyzed using 
the Ackland- Jones parameter [16j , which uses bond angle 
distributions to classify local environments as BCC, FCC, 
or HCP. The accuracy of this analysis was enhanced by 
position-averaging over 800 MD steps. Precipitate struc- 
tures have been rendered using OVITO [17] . 

Figure [T] exhibits precipitate structures observed in 
simulations at 700 K. Initially, Cu precipitates are 
isostructural with the surrounding BCC-Fe matrix, as 
shown exemplarily in Figure [l|a) by a precipitate con- 
taining approximately 15,000 atoms. Even larger pre- 
cipitates undergo a structural transformation from the 
BCC phase to a multiply-twinned 9R structure as il- 
lustrated in Fig. [ljb) where a cluster containing about 
36,000 Cu atoms is shown. The 9R phase is a closed- 
packed (CP) lattice that differs from FCC through its 
stacking sequence of the CP planes, as shown on the 
right-hand side of Fig.[ljc). Using the Ackland- Jones pa- 
rameter for structure identification (see Method section) , 
we obtain sequences of FCC and hexagonal close-packed 
(HCP) atoms in Fig. [ljb) that are in accord with the 
ideal 9R structure, c.f. Fig. [l](c). In addition, we observe 
stacking faults as well as several twins as indicated by ar- 
rows in Fig.jljb). These results closely resemble the TEM 
micrographs of 9R Cu precipitates in cv-Fe [U EJ 0] . Fur- 
ther analysis reveals that also parameters such as twin 
spacing and misorientation angles agree very well with 
experiments. 

In order to map out the size-temperature phase dia- 
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FIG. 1. (a,b) Cross sectional (top) and perspective (bottom) views of (a) BCC and (b) twinned 9R Cu precipitates at 700 K. 
The precipitates contain approximately 15,000 and 36,000 atoms, respectively. The 9R precipitate in (b) exhibits a pronounced 
BCC wetting layer, (c) Schematic of 9R structure [HE] which appears as a periodic sequence of two HCP atoms and one FCC 
atom along the [001]gR direction when using the Ackland- Jones parameter. 



gram of Cu nanoprecipitates, we record the number of 
Cu atoms in CP/BCC local environments as a function 
of precipitate size at various temperatures. The open 
symbols in Fig. [2]ja) show the evolution of these popu- 
lations with increasing precipitate size at 700 K. Precipi- 
tates exceeding 28,000 atoms are found to undergo a dis- 
continuous structural transition as they transform from 
BCC to multiply-twinned 9R structures. To study the re- 
verse transition, we start from a 9R precipitate obtained 
from an earlier simulation and gradually reduce the Cu 
content in the simulation. The thus obtained data are 
shown by the filled symbols in Fig. [2|a) demonstrating 
that at approximately 18,000 atoms precipitates revert to 
the BCC structure. There is thus a pronounced hysteresis 
between forward (BCC^9R) and backward (9R^BCC) 
transitions at 700 K extending from about 18,000 atoms 
to 28,000 atoms. 

The above analysis was repeated for six temperatures 
between 200 and 700 K, the results of which lead to the 
size-temperature phase diagram depicted in Fig. [3] It 
demonstrates that the critical cluster size for the BCC- 
9R transition is strongly temperature dependent increas- 
ing from approximately 2,700 atoms (3.9 nm diameter) at 
200 K to about 22,000 atoms (7.9 nm) at 700 K. The tran- 
sition is decidedly first-order at high temperatures with a 
pronounced hysteresis that vanishes below 300 K, where- 
upon the transition becomes nearly continuous. This 
is depicted in Fig. [2jh), where the numbers of BCC 
and CP atoms change around the critical size of 2,700 
atoms reversibly with negligible discontinuity. A more 
detailed inspection of the precipitate structure at 200 K 
(see movies and snapshots in Supplementary Material) 
shows that precipitates below the transition point (2,700 
atoms) have completely adopted a BCC lattice structure 
while precipitates containing more than 3,000 atoms ex- 
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FIG. 2. Number of atoms in different local environments at 
(a) 700 K and (b) 200 K. Filled and open symbols in (a) show 
results from forward (BCC— s>9R) and backward (9R— >BCC) 
simulations, respectively. The BCC-9R transition exhibits a 
hysteresis at higher temperatures that vanishes as tempera- 
ture is lowered. 



hibit multiply-twinned 9R structures as described above. 
In the intermediate size range [region A in Fig. |J»] the 
atoms that constitute the core of the precipitate rapidly 
fluctuate between different CP motifs. These fluctuations 
are absent at higher temperatures where the BCC-9R 
transition is sharp. 
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FIG. 3. Diagram representing the structure of Cu precipi- 
tates as a function of temperature and size. The light colored 
regions next to the phase transition line delineate the extent 
of the hysteresis. The hysteretic critical point is located along 
the dashed line between 200 and 300 K. 




FIG. 4. (a) Potential energy of atoms in Cu precipitates as 
a function of precipitate size for different temperatures. The 
vertical black bars represent the latent heat AUtr associated 
with the transformation, which is plotted in (b) as a function 
of temperature alongside the entropy change at the transition 
(AStr = AUtr/Ttr)- The upward (downward) pointing tri- 
angles show the number of core (wetting layer) atoms at the 
transition. Lines are guides to the eye. 



It is important to note that irrespective of tempera- 
ture even after the transition to 9R the precipitates still 
contain a substantial number of BCC atoms, which grows 
with precipitate size, albeit much slower than the number 
of CP atoms. This observation reflects the existence of a 
BCC-Cu wetting layer between the BCC-Fe matrix and 
the CP core of the Cu precipitate. It can be as thick as 
4 atomic layers at the transition point at 200 K, while it 
only amounts to about 1 to 2 atomic layers for very large 
precipitates. The BCC-Cu wetting layer is also clearly 
visible in Fig. [ljb) . In the following, we will see that the 
wetting layer plays a crucial role in the BCC-9R transi- 
tion, particularly at low temperatures. 



An important signature of first-order transitions is the 
release of latent heat, which is absent in continuous tran- 
sitions. To conclusively resolve the character of the ob- 
served transitions we have extracted the internal energy 
density of the Cu atoms in the precipitates as a func- 
tion of precipitate size at various temperatures as shown 
in Fig. [4ja). At high temperatures the data exhibit 
pronounced discontinuities when the higher-energy BCC 
phase transforms to the more stable 9R phase. As the 
temperature is lowered the discontinuity in the poten- 
tial energy is reduced until it apparently vanishes be- 
tween 200 and 300 K, see Fig. |4|b). If true, this implies 
that a critical point exists in this temperature range, be- 
low which the BCC-9R transition occurs continuously. 
However, the observation of a continuous transition con- 
tradicts bulk thermodynamics, according to which the 
distinct lattice symmetries require the BCC-9R trans- 
formation to be first order. Can this rule be violated for 
precipitate sizes as astoundingly large as 2,500 to 3,000 
atoms? 

The key to answering this question lies in the het- 
erogeneous distribution of energy and entropy inside Cu 
precipitates. For the purpose of this discussion we es- 
timate entropy by employing the Debye model for the 
vibrational spectrum, which provides a simple relation 
between entropy and mean square displacement (MSD). 
[IB] Figure [5]ja,b) shows the average atomic potential en- 
ergies (PE), and Fig. [5^c,d) depicts the MSD as a func- 
tion of distance to the centers of precipitates at 700 K 
for several BCC and 9R precipitates. The PE and MSD 
profiles of the BCC precipitates have similar shapes with 
a constant region in the core that remains unchanged as 
the precipitates grow. The wetting layer is identified as 
the region outside the core where both the PE and the 
MSD of the Cu atoms drop steeply. This confirms that 
the vibrational entropy of the core region of BCC precip- 
itates is much larger than the wetting layer entropy. It 
leads to the stabilization of BCC precipitates at elevated 
temperatures and derives from the vibrational modes of 
bulk BCC-Cu with imaginary frequencies that are stabi- 
lized through confinement by the Fc matrix. 

Inspection of the PE and MSD profiles of the 9R pre- 
cipitates in Fig. [5jb,d) reveals a more complex behavior 
than in the case of the BCC precipitates. While the wet- 
ting layer seems to be quite similar in size, energy and 
entropy to the BCC precipitates, the core now consists 
of two regions: (i) the inner core with uniform energy 
and entropy, and (ii) the outer core, which constitutes 
the interface region between the twinned 9R inner core 
and the wetting layer. It is strained to accommodate the 
misorientation between the inner-core 9R lattice and the 
BCC-Fe matrix. In this region, both energy and MSD 
increase but never exceed the values of the core of the 
BCC precipitates. Hence the 9R core always contributes 
positively to the latent heat. However, the total latent 
heat AUtr generated by the structural transformation is 
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FIG. 5. (a,b) Potential energy and (c,d) mean square displacement as a function of distance from precipitate center for (a,c) 
BCC and (b,d) 9R precipitates at 700 K. (e) Latent heat AUtr [from Fig. |4|b)] as well as core At/core and wetting layer j w i 
energy differences between BCC and 9R precipitates obtained by analyzing PE profiles. 



always smaller than the core contribution since the wet- 
ting layer energy is always slightly lower in BCC than in 
9R precipitates. The total latent heat and its associated 
entropy AS tr = AU tr /T are shown in Fig. fflh). Both 
decrease dramatically upon cooling mainly because at 
lower temperatures transitions involve much fewer atoms 
as they occur at much smaller cluster sizes. Furthermore 
at lower temperatures the latent heat per atom drops sig- 
nificantly [see Fig. [5je)] due to the fact that in smaller 
precipitates the core region, which provides the positive 
contribution to the latent heat, shrinks while the com- 
pensating wetting layer thickens. This is supported by 
the number of atoms in the core and the wetting layer 
regions at the transition as depicted in Fig. [4jb) . 

Nevertheless based on the reported temperature de- 
pendence of the entropy of transformation in Fig. Qb) , 
one expects the BCC-9R transformation to involve a di- 
minishing but finite latent heat at temperatures clearly 
below the critical temperature (« 300 K) at which the 
hysteresis vanishes. Note that at this temperature, as 
many as 500 Cu atoms form the core of the 9R precipitate 
and thus directly participate in the transformation, see 
Fig. [4jb). These observations thus support the existence 
of a dynamic critical point in the BCC-9R transition 
of Cu nanoprecipitates. At this temperature, the nucle- 
ation barrier between the two phases disappears. The 
microscopic origin of this behavior lies in the mechani- 
cal instability of bulk BCC-Cu on the one hand and the 
stabilization of a BCC wetting layer with stiff phonon 
modes at the interface with the BCC-Fe host lattice on 
the other. 

It is important to note that while our observations can- 
not prove the existence of a thermodynamic critical point 
(TCP) in the BCC-9R transition of Cu nanoprecipitates, 
they suggest that the phase diagrams of nanoprecipitates 
may contain TCPs that have no counterpart in the corre- 
sponding bulk thermodynamic limit. It is of course well- 
known that the phase diagrams of nanoparticular systems 



deviate from the bulk due to confinement. In particu- 
lar the loci of phase boundaries become size dependent 
and thus phase transition temperatures and pressures are 
modified for liquid-gas |19| . liquid-solid p2PH!24"] as well 
as solid-solid transformations 25-34 . Furthermore, dy- 
namic critical points leading to vanishing hysteresis have 
been observed to precede the TCP in liquid-gas transi- 
tions [19j [35] . In contrast, one might expect the first- 
order character of transitions between phases with dis- 
tinct symmetries to be conserved. The presence of a TCP 
in BCC-9R transformations of Cu nanoprecipitates vio- 
lates this expectation. It is caused by the heterogeneous 
nature of nanoprecipitates, where the total latent heat 
can vanish with the core and wetting layer regions ex- 
hibiting compensating discontinuities. This heterogene- 
ity is best illustrated in Fig. [5^e), which shows the total 
latent heat per atom, as well as the inner core and wet- 
ting layer energy differences between BCC and 9R pre- 
cipitates. For temperatures above 500 K, the contribu- 
tion to the latent heat from the inner core region reaches 
its asymptotic limit. Below this temperature, this con- 
tribution is slowly reduced as the inner core region is 
gradually shrinking and becoming strained. Figure [5je) 
demonstrates that the latent heat per atom can nearly 
vanish while significant differences in the energetics of 
BCC and 9R precipitates remain. This points to the fact 
that the balance between wetting layer and core ener- 
gies is fundamental to the low-temperature criticality in 
this system. Furthermore, if a TCP existed, higher order 
derivatives of the free energy would become discontinu- 
ous rather than diverge. 

The phenomenology described in this work not only 
has direct consequences for our understanding of the Fe- 
Cu system but transpires to many more materials. Based 
on the analysis presented here one can anticipate the oc- 
currence of similar effects in other immiscible systems, 
whenever the minority phase is unstable in the lattice 
structure of the host. 
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